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Simple Monte Carlo Computer Procedure for the Depth Parameters 
Determination of Implanted Ions in Amphous Targets 


Prosta metoda Monte-Carlo okreslania parametrow rozkladu implatowanych jonow 
w amorficznych tarczach 


Модифицированный метод Монте-Карло для определения параметров распределения 
имплантированных ионов в аморфные мишени 


INTRODUCT ION 


The Monte Carlo simulation of the slowing down of energe- 
tic ions in solids is now widely used in studies of ion im- 
plentation, radiation damage and sputtering. It allows more 
rigorous treatment of the depth profiles of implanted ions 
giving all required higher moments of the distribution after 
simple evaluation. This is rather difficult when the present 
analytical formulation based on transport theory is used. The 
major inconvenience of the Monte Carlo nethod is that it requi- 
res large computer-time usage, especially for high energy ions 
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treatment with suitable statistical precision. Each simplifi- 
cation of the method reducing the computer time consumption is 
then important. 

There are several ion transport procedures based on the 
computer simulation of the slowing down of energetic ions in 
solids [1 - 11] . One of their major differences lies in the 
treatment of the nuclear scattering. For example, Robinson et 
al. [2 - 4] treat this scattering by numerically evaluating 
classical scattering integral. Other authors used the momentum 
approximation extended to large angles  [5, 6, 10] ог trunca- 
ted Coulomb potential [7 - 9] . Biersack and Haggmark [11] 
base their formalism on a simple but precise analytical ex- 
pression for determining nuclear scattering angles for Molier 
potential /TRIM - program/. The TRIM program reduces computer 
usage by at least эп order of magnitude in comparision with 
procedures based on the precise numerically treated classical 
scattering integrals. 

In the paper we propose the Monte Carlo procedure based 
on an alternative analytical method of evaluation of the nuc- 
lear scattering angles for the Molier potential. The program 
reduces significantly computer time consumption also in com- 
parision with fast TRIM procedure keeping up nearly the same 
eccurancy of evaluation, 


1. NUCLEAR ANO ELECTRONIC ENERGY LOSS EVALUATION 


During the slowing down of energetic ions in solids, the 
particle undergoes many collisions gradually losing the energy 
as а result of elastic /nuclear/ and inelastic /electronic/ 
interaction with target's atoms. The nuclear energy lose in a 
single collision can be easily calculated when scattering an- 
gle of the projectile is known. Last one can be treated from 
clasical scattering integral 
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where b is e impact parameter, E, is the relative kinetic 
energy, г і the interatomic separation, Vir) is the potential 
of interatomic force, and К is the apsis of the collision de- 
fined by [:- 9% 0 

0 Er / 

The scattering angle can be also evaluated from approxima- 
tion analytic formula for example proposed by Biersack and 
Haggmark [11] . 

In this program we use the analytic equation of scattering 
angle evaluated from Lindhard's [12] nuclear differential 
cross section of the form 


dT 
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where B is the Bohr's collision diameter,a is the screening 
length, kn is а constant, T the energy transfer in the collisions 
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M, and M, are the masses of the projectile and target atom res- 
pectively, Е the initial kinetic energy of the projectile and n 
is the parameter of the inverse power potential [13] . 

By the integration of equation (1) the energy transfer in 
the collision can be evaluated as 


b? Е А -n 
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where b is the impact parameter, and 
Е ‚1407 MM 
А, = 1,4397: 10 um 2122 lev. ст] 


On the other hend 


T = Tm sin? — 
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where Ka is a new constant 


Р = b پر‎ is the reduced impact parameter 
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is the reduced energy. 
In the calculation, we have chosen to use the Firsov 
screening length given by 


where AL 0,529 Я is the Bohr radius. 
In the equation /2/ parameter n » 1 corresponds to Ruther- 
ford scattering п = 2 to weak screening collision and n = 5 
to hard sphere collision. In the Monte Carlo simulation method 
equation /2/ cannot be directly used /with constant value of 
n and K, because each “history” of an ion begins with high e- 
nergy and particle undergoes many collisions before it stops. 
The п as a function of £ and к. as a function of P and € 
/from the fitting to the Molier potential (14] / can be repre- 
sented by the following expressions 


n(£) = 1*4 exp 19.8601) i (3) 
and 
Kn = A-D-P-C * 


where 


pe" \ 
c? )] 51:295 an ( ] 


A = 0,224 [ 1. 45210 (1+ | 13,8 +0.12 VE 


Simple Monte Carlo Computer Procedure .,, 21 


= * 
| (z/e),us © e 


a 
42 
сч 
©. 
сч 
€ 
| ч б 
| з 
| 5 
— = — A S 1 
| 5 5 8 
| (2/01 „мі o 


Fig. 1. 112 (9%) 88 а function of reduced impact parameter P 
for reduced energy between 10^? and 10 
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Figures 1 а, b show the comparison between sin (95) values 
calculated from the equations (3) and  (4)and from the exact 
evaluation of scattering integral for Molier potential. Formu- 
lation presented above simplifies evaluation of scattering angle 

8 in comparison with TRIM procedure while at the same time 
sacrifices a little accuracy. 

Inelastic energy losses are included in the program by 
using the model in which electronic energy loss is treated in- 
depedent of the nuclear energy loss. In this case, neglecting 
impact parameter dependence, relation between electronic energy 
loss and distanca L travelled between collisions can be герге- 
sented by 

Te = LNSẹlE) 


where S /Е/ is the electronic stopping cross section, To cal- 
culate $е/Е/ we use Lindhard ~ Scharff [15] model 


Se lE! = ке 2 
where 


1212 216 72 
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Оп the other hand Оеп and Robinson [16] have suggested, 
that the electronic energy loss as a function of closest 
approach can be expressed as 


Te = 045k VE .- 03 Ro/a 


Tae (6) 
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To save computer time usage we introduce also approximation 
formula for Во in the form 


Yn 
0,15 
[Ro lapp*a 22 s b (7) 


Application of the approximation /7/ to equation /6/ reduces 
the inelastic energy losses below the prediction of Lindhard 
[15] at low energies and small impact parameters but not to 
such and extent as predicted by Oen and Robinson [16]. This is 


easily seen in the Fig. 2 which compares the electronic stoppinc 
cross section 
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Fig. 2. Electronic stopping cross sections for T incident on Cv 
in different approximations. Curve /с/ - according to equation 
/6/ with the approximation /7/ 
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for T incident on Cu calculated according to the Lindhard-Scharff 
and Oen-Robinson models end according to equation (6) with the 
approximation 17} . Equation 16] has been made available as 
an option in our program. 


2, DISTANCE BETWEEN COLLISIONS 


The targets were assumed to be amorphous, so that the dis- 
tance between collisions L can be assumed as equal to the mean 
atomic separation L s м-1/з where N is the stomic density of the 
target. This assumption leads /according to Biersack and Haggmerk 

[11] / to the determination of impact parameters 


е бшш 
b = [R Z N 72 


where R, is the random numbers which are evenly distributed 
between О and 1. The value of L = N71/3 is not Strictly valid 
at low enorgies where the distance between collisions depends 
upon the impact parameter b and mass ratio м = М,/'1;, Іп general 
mean free path becomes shorter than the distance lto the next 
target atom since the deflection point of projectile atom tra- 
jectory lies in front of the scattering centre. This path - 
length reduction /indicated as X4 in the Fig. 3/ corresponds 
to the so called “time integral T " which is given by the ex- 
pression similar to the scattering integral [4]. 

The value of x, can be expressed by 


x, 2 E -') d 6 71, U. 7 


In the low energy region, assuming the hard sphere collisions 
model, the time integral T can be approximated by its hard-core 
value, T = btg 8/2. Such assumption leads to 


xz b tg V 


and 
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The procedure presented above has been used by Biersack and 


Haggnark (11] in the TRIM program and is also adopted in this 
work. 


INITIAL TARGET 
ATOM LOCATION 


Fig. 3. Trajectories of two interacting particles 


At high energies only few of the many collisions cause signi- 
ficant deflection from the straight path of flight. The program 
is set up to eelect the smallest impact parameters and neglect 
the other collisions with larger impact parameters for which 
sin? 8/2 becomes less then — 10 4. We find that the smallest 
reduced impact parameter can be approximated by the formula 


"nin (253-1) for £400 


For the parameters larger than Pane scattering angles of pro- 
jectile 8 and nuclear energy loss T, are assumed to O and 
only electronic energy loss is taken into account. 
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3. RESULTS AND COH PARIS IONS 


To test the extremely simplified computational technique 
presented here we use the same set of experimental and theo- 
retical data as proposed by Biersack and Haggmark [11] . In 
the Fig. 4 our results for mean projected range В and stan- 
dard deviation 6 for 20-180 keV boron implanted into silicon 
is shown together with the experimental data of Hofker et al 
17 Jend Ryssel, et al [18-] and with theoretical results of 
Biersack and Haggmark /TRIM program 11) / and Brice £190. 


THIS. WORK 


Р . rere: 
50 100 150 200 
ENERGY [kev] 


Fig. 4. Mean projected range R_ and standard deviation € for boron 


implanted into silicon /0 - experimental data of Hofker at а1.[17], 
О - experimental data of Ryssel et al, [18] , — = Brice [19] 


- k/ky = 1.59,x- TRIM - к/к = 1.59, А = this work k/kL = 1.59/ 
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The follwing definitions have been used for moments calculations 
from the Monte Carlo histograms 


<%= Rp = jx fix) dx 


ftti dx 


1 
а fix) х e 
Ji dx 


/mean projected range/ 


/ depth straggle/ 


Ig AE fix) dx 


/skewness/ 
tix) dx 


1 
59 


к Six -<x м“ f(x) dx 


E п тата /kurtosis/. 


In the Figures 5 and 6 projected ranges of antimony in 
silicon and helium in copper calculated by Biersack and Hagc- 
mark /TRIM/ and Oen and Robinson [20] - /MARLOWE program [4] / 
are compared with results obtained in this work. 


PROJECTED RANGE [А] 


ENERGY [keV] 


Fig. 5. Projected ranges of Sb ions in silicon /о = experimental 
data od Oetzman et al. [21] , — = TRIM results LIST „ A 
this work/ 
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Fig. 6. Mean projected ranges of low energy helium ions implan- 
ted into Cu / ~= = MARLOWE , ® = TRIM, x this work (T. = 


according to Oen and Robinson [16] with Ry approximation 
/equation 7/)/ 


In the case of antimony ions in silicon the experimental data o! 
of Oetzman [21] аге included in the figure for comparison. 

In the Fig. 7 the Monte-Carlo reduced sean projected depths 
estimated by Latta [10] for Bi implanted into Ge are compared 
with present calculation and some of the experimental data [21, 
22, 23]. In this case the mean projected ranges are represen- 
ted in reduced units by the formula 


Ra · 4 Nae 2 
eie ee МММ ама 


where 


-% 
атр= 0.8853 · 0028 + 2,3) 


and а = 0,529 À 
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Fig. 7. Reduced mean projected range фр versus reduced energy 
Е, /Experimental and theoretical data compilation of Latta 
[ 10]/Experimental data: O- Oetzmann et al. [21] ‚ А Grant 
et al. [22] , o - Besenbacher et al. [23] . 


We also compare the values of the higher moments /skewnes and 
kurtosis/ with those calculated by Latta [10] . The results 
are summarised in Tab. 1. 

The results of the calculations presented above indicate 
rather good agreement with other theoretical results and also 
with experimental measurements except the heavy ions implanta- 
vion /see Fig. 5 and 7/. In this case the theoretical predic- 
tions are lower than those of the experiment, All theoretical 
curves /Fig. 7/ converge at the largest Е value. At low € 
values the best fit to the experimental data are obtained when 
Monte-Carlo treatment for Solid State Thomas-Fermi /SSTF/ 
potential /Latta [10] / is used. In general, neither Thomas- 
-Fermi potential nor SSTF and Molier approximation to the 
Thomas-Fermi screening function, results in mean-depth estimates 
that would agree with experiment over the full € range [10, 11] . 
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Tab. 1. Comparison of Monte Carlo yt / 
results with Minterbon's [25] numerical estimates at reduced 
energy of Е = 0,01 


Collision MC-Latta MC-This Numerical 


y bie ^ pair йо] work Winterbon [25] 
Mean depth В1=> Ge 102 84,5 67,6 
«x» /А/ Ars» Au 18,7 16,2 17,1 
Depth straggle Bis» Ge 38,5 35,2 25,7 

8 „R Аг= Au 2807 19,5'V 19,2 
Skewness Biz» бе 0,485 0,539 0,546 

S Ars? Au O20 — ОШ 0,23 
Kurtosis Biz» Ge 3,32 3,43 3,23 

K Aras Au 3,21 3,487 3,18 
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Xx - calculated for infinite medium 


4. CONCLUSION 


In the paper we have presented the simple Monte Carlo simu- 
lation procedure of slowing down of energetic lons in emorphous 
materials, One of the features of this program which distin - 
quishes it from other procedures is the use of extremely simpli- 
fied analytic expression of scattering angle and energy loss 
evaluation. The other one is connected with electronic stopping 
cross section treatment, in which we propose an alternative 
approximation to the impact parameter dependence of the elec- 
ctronic energy loss. 

To save computer time usage in the cese of high energies 
calculations the simple impact parameter selection procedure 
is presented. 
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The agreement between results presented here and those 
based on TRIM program is rather sotisfactory. The computer time 
consumption, in the comparison with TRIM procedure is lowlering 
significantly. 
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STRESZCZENIE 


W pracy przedstawiono uproszczona методе Monte-Carlo okre- 
$lenia parametrów gleboko$ciowych rozkładów jonéw implantowa- 
nych do tarcz amorficznych. W obliczeniach przeprowadzonych 
dla kilku kombinacji jon-tarcza wykorzystano potencjai Moliera. 
Otrzymane wyniki porównano z obliczeniami teoretycznymi i wy- 
nikami doświadczalnymi innych autorów. 


РЕЗЮМЕ 


В работе представлен модифицированный метод Монте-Карло 
определения концентрационных распределений имплантированных 
ионов в аморфные мишени. Расчеты проводились для нескольких 
пар ион-мишень с потенциалом Мольера. Полученные результаты 
сравниваются с теоретическими предсказаниями и эксперименталь- 


ными данными других авторов. 
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